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SUMMARY 


Various  aspects  of  ducted  propeller  theory  and  design  are  considered. 
These  include  the  lift  of  a  moderately  inclined  ducted  propeller, 
possibilities  for  increasing  the  static  thrust,  the  related  blading 
design  and  flow  field  analysis  of  interference  with  adjoining  wings 
or  bodies . 

The  inclined  duct  theory  is  shown  in  agreement  with  several  sets  of 
test  data  from  different  sources. 

Large  static  thrust/horsepower  values  with  low  jet  velocities  are 
predicted  for  designs  which  accelerate  the  inflow  in  accordance 
with  specified  pressure  changes  inside  the  duct.  These  consid¬ 
erations  are  at  present  without  experimental  verification.  For  such 
flows,  a  range  of  appropriate  blading  designs  are  shown  by  means  of 
solidity,  pitch  distribution,  jet  velocity  and  tip  speed. 

Digital  computer  studies  are  recommended  to  evaluate  favorable 
and  unfavorable  interference  flow  arrangements  between  ducted 
systems  and  surrounding  surfaces.  The  fan-in-wi^g  flow  field  with 
a  jet  of  finite  size  is  one  of  many  cases  which  can  be  handled. 
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LIFT  ON  AN  INCLINED  DUCTED  PROPELLER 


A.  General 

It  is  well  known  that  a  duct  of  length  L  and  inner  diameter  d  inclined 
at  angle  a  to  an  oncoming  flow  produces  a  lift  force.  The  lifting 
characteristic  of  the  inclined  duct  (ring-wing)  is  similar  to  that  of 
the  conventional  planar  wing;  this  is  seen  from  the  aspect  ratio  type 
of  effect  dependent  on  L/d  and  the  lift-drag  polar.* 

Analytical  studies  now  enable  prediction  of  the  lift  curve  and  the  de¬ 
pendency  on  L/d  to  a  degree  quite  useful  for  engineering  prediction.** 
This  is  demonstrated  in  REF.  1.  It  is  worth  noting  that  the  use  of 
the  inner  diameter  of  the  duct,  for  ducts  of  finite  thickness  ratios, 
seems  to  provide  the  favorable  agreement  between  theory  and  ex¬ 
periment.  These  results  are  established  for  ducts  of  essentially 
constant  inner  diameter.  Similar  developments  for  ducts  having 
appreciable  inner  contour  change  are  not  clearly  evident’ in  the  litera¬ 
ture. 

With  the  introduction  of  an  operating  propeller  inside  the  duct  the  lift¬ 
ing  characteristics  become  different  from  the  unpowered  duct.  From 
an  elementary  viewpoint  this  is  seen  from  the  larger  vertical  flow  (jet) 


*  REF.  1  Experimental  Investigation  of  Lift,  Drag  and  Pitching 
Moment  of  Five  Annular  Airfoils  by  H.S.  Fletcher, 
NACATN  4117.  1957. 


**  REF.  2 
REF.  3 
REF.  4 


The  Ring  Airfoil  in  Nonaxial  Flow,  by  H.S.  Ribner 
Jour.  Aero.  Sci. ,  Vol.  14,  1947,  page  529. 

Etude  Theorique  De  L'Aile  Annulaire  by  M.G.  Faure. 
Technique  Et.  Science  Aeronautique.  Tome  6,  1956. 

Zur  Aerodynamic  des  Ringflugels  by  J.  Weissinger, 
DVL  Bericht  Nr.  2,  1955. 
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momentum  component  of  the  powered  versus  the  unpowered  duct. 

REF.  5,  for  example,  uses  this  type  of  approach  together  with  empir¬ 
ical  factors  to  develop  a  method  to  predict  the  lift  force. 

A  lift  dependence  on  L/d  for  the  powered  duct  can  also  be  anticipated. 

The  extent  that  the  lift  departs  from  the  vertical  component  of  jet 
momentum  depends  on  L/d.  For  example,  a  jet  engine  enclosed  in  a 
nacelle  will  have  an  L/d  value  in  the  order  of  3;  inclining  the  engine, 
at  relatively  low  airspeeds,  yields  a  lift  force  very  close  to  the  ver¬ 
tical  jet  exit  momentum  component.  With  the  introduction  of  the  ducted- 
fan  jet  engine,  the  enclosing  nacelle  shows  smaller  L/d  values  and 
it  can  be  anticipated  that  the  interaction  of  internal  and  exit  flows  will  pro¬ 
duce  vertical  forces  greater  than  the  jet  momentum  component.  An 
introductory  treatment  of  the  effect  of  duct  length/diameter  ratio  in 
relation  to  interference  effects  between  duct  and  jet  flow  is  given  in 
REF.  6.  Increasingly  large  L/d  values  lead  to  smaller  interaction 
effects. 

In  the  next  section,  available  test  data  on  moderately  inclined  powered 
ducts  are  compared  with  an  analytical  procedure  based  on  flow  through 
nonpowered  ducts. 


REF.  5  Performance  Calculation  of  Aerodyne  Systems  in 
Cruising  Flight  by  A.M.  Lippisch.  Collins  Radio 
Company  Rept.  Cer-617,  1957.  Confidential 

REF.  6  An  Estimate  of  the  Forces  on  Annular  Fairings  by 

J.A.  Bagley.  Jour.  Roy.  Aero.  Soc.,  vol.  63,  1959, 
page  315. 
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B.  Lift  of  a  Moderately  Inclined  Ducted  Propeller 

Details  of  the  method  are  described  in  Appendix  A.  This  method,  as 
shown,  extends  the  vortex  theory  of  an  unpowered  duct  to  take  account 
of  the  induced  effects  of  the  jet  vorticity. 


For  the  unpowered  duct  the  lift  coefficient  ^  is  given  by 


C.  (a0/ 5  7.3) 


(1) 


where  is  the  lift  curve  slope  dependent  on  the  duct  length  diam- 

Ct 

eter  ratio.  Referenced  papers  present  different  methods  for  deriving 
the  lift  curve  slope.  For  simplicity  the  curve  given  in  REF.  1  can  be 
taken  in  the  approximating  form: 

CL  =  [.0547/(L/d)-.0063]  /(L/d)2  .  (2) 

a,u 


Using  this  relation  for  the  lift  curve  slope,  the  derivation  of  Appendix 
A  shows  the  lift  coefficien 
ducted  propeller:  namely: 


A  shows  the  lift  coefficient,  ^  for  the  duct  with  an  operating 


C  =  {[.25Cl  {\N+  l)2]  +(.57T/\)[(V./V)(V1/V)-l]}a°/57.3  (3) 

’  P  m  J 


whe  re 


V  is  the  averaged  velocity  inside  the  ducted  system. 

V  is  the  advance  speed. 

V.  is  the  jet  velocity. 

X  i  s  L  /d  • 


Without  power  ,  the  theory  requires  V  =  V  and  V.  =  V,  and  Eq.  3 

J 

r  educed  to  Eq.  1  for  these  velocity  values. 


-3- 


An  indication  of  the  effect  of  power  on  the  lift  curve  slope  can  be  had 

by  setting  V  /V  =  1.0  in  Eq.  3  and  inserting  Eq.  2  for  the  unpowered 

CjJ  ;  the  result  becomes 

a, 


_°LlR 


a,u 


1  +  (.57 T/\)  (V./V-1)/Cl 
J 


(4) 


and  is  plotted  on  FIG.  1  with  the  thrust  coefficient  C^,  as  abscissa; 
the  definition  of  C^,  is 


CT  =  T  / ( 1  /2)pV2  Aex 


(5) 


where  A  is  the  exit  area  of  the  duct.  For  the  assumption  that  the 
ex  r 

ultimate  jet  area  A.  =  A  ,  it  follows  that  for  a  radially  averaged  flow 

J 


CT  =  2(V/V)(V ,/V-l)  (6) 

J  J 

based  on  the  thrust 

T  =  pA.V.(V.-V),  (7) 

J  J  J 

FIG.  1  shows  a  large  increase  for  the  lift  slope  curve  with  increasing 

thrust  coefficient.  Decreasing  L/d  ratio  also  increases  the  lift 

curve  slope.  For  L/d  values  of  0.6  and  greater  the  effect  of  L/d 

becomes  much  smaller.  The  valid  extent  of  FIG.  1  depends  on  the 

closeness  of  A  to  A.,  With  increasing  thrust  coefficients,  it  is 

e  j  6 

most  likely  the  jet  contracts  downstream  of  the  exit  plane.  Jet  con¬ 
traction  probably  is  a  secondary  effect  for  thrust  coefficients  of  1.0 
and  below. 


FIG.  1  is  based  on  V  /  V  =  1.0;  larger  lift  slope  curves  will  occur 
for  designs  having  V1/V>1.0. 
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CL,p  1,U 

13.' 


11.  Ol 


CT=T/(i)pV2Aex 
-  Lift/(l/2)pV2  Ld 


Thrust  Coefficient,  C, 


FIG.  1  RATIO  OF  LIFT  CURVE  SLOPES  FOR  POWERED 
AND  UNPOWERED  DUCTS 
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C.  Comparison  with  Test  Data 


Eqs.  1  and  3  will  be  compared  with  data  in  REFS.  7,  8  and  9.*  In 
order  to  compare  theory  with  experiment,  it  is  necessary  to  know 
flow  velocities  at  the  duct.  Such  data  are  not  available  for  the  refer¬ 
enced  experiments  and  the  procedure  outlined  below  is  followed  for 
estimating  the  jet  velocity. 


The  power  input  to  the  flow  (ft.  -lb.  /sec. )  according  to  kinetic  energy 
principles  is; 


P  =  (l/2)pA.V-’(V/V)[(V./V)2  -l]k  (8) 

J  J  J 


where  V 


V. 

J 

A. 

J 


is  the  free  stream  velocity. 

is  the  jet  velocity. 

is  the  jet  cross  section  area. 

is  a  loss  function  greater  than  1.0;  for  zero  loss,  k^  = 


1.0. 


The  propulsive  efficiency  of  the  ducted  system  from  momentum 
consideration  is: 


e  =  2/[kT(V/V+  1)] 
P  L  J 


(9) 


where,  as  above,  k^  takes  account  of  losses  through  the  blading. 


*  REF.  7  Shrouded  Propeller  Investigations:  Wind  Tunnel  Tests 
of  a  Shrouded  Propeller  with  a  10-Bladed  Propeller, 
Exit  Stators  and  Long  Chord  Shroud  With  High  Speed 
Inlet  and  No  Exit  Diffusion  by  V.O.  Hoehne.  Univ.  of 
Witchita,  Engrg.  Rept.  No.  213-5,  Jan.  1959. 

REF.  8  Wind  Tunnel  Test  of  Several  Ducted  Propellers  In  Non- 
Axial  Flow  by  W.  J.  Gill,  Hiller  Aircraft  Rept.  No. 
ARD-224,  April  1957. 


REF.  9  Wind  Tunnel  Tests  of  Shrouded  Propellers  At  Mach 
Numbers  Fropi  0  to  0.60  by  R.  M.  Grose.  WADC 
TR  58-604.  December  1958. 
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Eliminating  kj  between  Eqs.  8  and  9  yields: 


<V./V)[(V/V)-1]  =  Pe  /pA  V3 

J  J  a  J 


The  right  side  of  Eq.  10  containe  P  and  e^  which  are  obtainable  from 

the  test  data;  jet  area  A.  can  be  approximated  by  the  duct  exit  area, 

J 

provided  the  jet  does  not  contract  greatly.  The  effect  of  jet  con¬ 
traction  is  considered  later. 


comparison  with  the  experiment  data  (REF.  7] 


REF.  7  presents  thrust  and  power  data  in  terms  of  coefficients 
defined  by 


CN  " 

F  /pV2  d2 
n 

T 

r 

T/pV2  d2 

(11) 

n*. 

ii 

P/pV3  d2 

V 

j 


FIG.  2  GEOMETRY  OF  FORCES 
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It  is  seen  from  FIG.  2  that 


CT  =  C-.cos  a  +  T  sin  a.  (12) 

L  N  c 

The  geometry  of  the  ducted  system  in  REF.  7  is  as  follows: 

a)  propeller  diameter,  13.92  inches. 

b)  hub  diameter,  4.76  inches  (d,  /d  =  .3). 

n 

c)  exit  diameter,  13.92  inches. 

d)  duct  length,  9.526  inches. 

e)  L/d  =  .  684. 

f)  exit  area  is  reduced  5.7%  due  to  support  fairing. 
TABLE  I  lists  the  test  data  to  which  Eq.  3  will  be  compared. 


TABLE  I 

Test  Data  From  REF.  7 
(light  loading) 


o 

a 

J 

o 

(3  @  ,8R 

P 

c 

T 

c 

e 

P 

Closer 

T  sin<* 
c 

L,p 

0 

1.35 

42.5 

.5 

.34 

.678 

0 

0 

0 

0 

10 

1.35 

42.5 

.5 

.34 

.678 

.3 

.296 

.059 

.  355 

20 

1.35 

42.5 

.5 

.34 

.678 

.6 

.564 

.116 

.  680 

30 

1.35 

42.5 

.5 

.34 

.678 

.8 

.  693 

.  170 

.863 

a  =  o° 

Inserting  the  geometrical  and  experimental  data  into  Eq.  10  yields 

V./V  =  1.36 
J 

Vj/V  =  1.28  (from  flow  continuity).  (13) 

These  values,  together  with  Eq.  3  yield  the  theoretically  determined 
value  of  Cj^  .  FIG.  3  compares  theory  and  experiment.  Also  shown 
on  FIG.  3  is  the  lift  curve  for  the  unpowered  duct. 
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Theory  is  shown  on  FIG.  3  to  furnish  a  correct  guide  to  the  physical 
situation  of  this  isolated  set  of  test  data;  for  moderate  angle  of  attacks 
the  theory  furnishes  a  good  engineering  approximation  despite  the 
simplfying  assumptions  involved.  The  addition  of  power  appears  to 
increase  lift  about  100%  within  the  unstalled  range;  this  ratio  may 
persist  even  to  the  higher  angle  of  attack  range  provided  the  lift  curve 
slope  of  the  unpowered  duct  decreases  similarly  to  that  of  the  powered 
duct. 


iLift 


Angle  of  Attack,  a,  degrees 


FIG.  3  LIFT  OF  INCLINED  DUCT  (L/d  =  .684) 
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2.  Comparison  with  REF.  7  including  the  effect  of  jet  contraction 


The  foregoing  comparison  of  lift  coefficients  is  made  for  a  ducted 
system  with  a  relatively  low  thrust  loading.  This  comparison  re¬ 
quires  an  estimate  of  the  jet  cross-section  area;  in  view  of  the  light¬ 
loading  conditions,  a  jet  area  the  same  as  ducted  propeller  exit  area 
is  assumed.  While  this  choice  is  reasonable  for  a  lightly  loaded 
system  it  should  not  be  expected  to  apply  generally  for  a  highly  load¬ 
ed  ducted  propeller  having  an  exit  area  the  same  as  the  propeller 
disk  area.  The  test  data  of  REF.  7  are  obtained  from  a  ducted  system 
with  only  a  slight  duct  contraction  between  propeller  disk  and  exit 
plane  due  to  support  struts.  Hence,  for  much  higher  thrust  loadings, 
the  flow  leaving  the  duct  must  be  presumed  to  return  the  equilibrium 
downstream  of  the  duct  by  an  externally  contracting  jet.  In  order  to 
compare  theory  and  experiment  for  duct  lift  for  higher  thrust  loadings 
at  low  propeller  advance  ratios,  an  estimate  is  needed  of  the  jet  area. 

To  illustrate  further  the  nature  of  this  problem,  consider  the  con¬ 
tracting  jet  at  the  duct  exit  as  shown  on  FIG.  4. 


FIG.  4  CONTRACTING  JET 

Bernoulli's  equation  between  stations  e  and  j,  together  with  the 
corresponding  continuity  relation,  yields,  for  an  incompressible  fluid, 

Ae/A.  =  V1  +  (pe’pj)/,1/2)pVe'  (14) 
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which  allows  an  estimate  of  the  contracted  jet  area,  provided  the 
overpressure  is  known  at  the  duct  exit;  this  pressure  can  be  pre¬ 
dicted  by  an  iterative  procedure  for  calculating  the  pressure  rise 
across  the  propeller  for  a  prescribed  power  input  and  advance  ratio. 
The  interacting  effect  of  the  duct  on  the  inflow  at  the  rotor  is  an  added 
complication.  One  general  approach  is  to  calculate  this  pressure  rise 
at  each  blade  radius  by  a  'strip'  method  taking  account  of  duct  inter¬ 
action.  From  the  pressure  rise  and  duct  cross-section  areas,  the 
pressure  at  the  duct  exit  can  be  estimated.  This  can  be  a  lengthy 
procedure  and  for  the  purpose  of  predicting  duct  lift  at  low  advance 
ratios  an  approximate  method  will  be  attempted  as  described  below. 

It  is  assumed  the  streamtube,  starting  from  the  duct  exit,  contracts 
similarly  to  the  jet  of  a  free  propeller  having  the  same  mass  flow  and 
ultimate  jet  velocity.  The  actuator  disk-momentum  method,  it  will 
be  recalled,  gives  the  velocity  through  the  disk  of  a  free  propeller  as 
the  average  of  upstream  and  downstream  velocities.  For  the  con¬ 
tracting  jet  leaving  the  duct  exit,  the  above  assumption  requires  the 
duct  exit  velocity  to  be  the  average  of  upstream  and  downstream  ve¬ 
locities.  That  is,  by  continuity  of  flow 

A./A  =  V  /V.  =  (1+V./V)/2(V./V)  (15) 

J  e  e  j  j  j 

which  furnishes  an  equation  for  the  additional  variable  A..  The  com¬ 
parison  between  theory  and  experiment  for  powered  duct  lift  can  now 
follow  the  method  applied  earlier.  A  convenient  modification  of  this 
method  can  be  based  on  calculating  the  jet  velocity  from  thrust  data 
given  by  the  coefficient  T£  (relations  11): 

T=  pV2d2Tc  =  pAV(V  -V)  (16) 

J  J  J 


-11- 


(17) 


together  with  Eq. 


15  for  A., 
J 


it  follows  that 


V./V  =  |/  1  +  2d2  T  /  A.  1 

J  |(  c  e 


from  which  the  jet  velocity  can  be  estimated  for  a  jet  contracting 
downstream  of  the  duct  exit. 


For  the  duct  geometry  of  REF.  7  shown  on  page  8,  Eq.  17  becomes 

V./V  =  n/1  +  2.967  T  .  1  (18) 

J  c 

In  addition,  the  internal  velocity  Vj/V  is  also  required  by  Eq.  3  for 
the  powered  lift.  From  the  continuity  of  flow  between  propeller  and 
exit  for  the  design  of  REF.  7,  the  following  ratio  applies: 

V./V  =  .943  (19) 

1  e 

and  from  the  exit  station  to  the  ultimate  jet 

v./v  =  ,943(A  /A  KV./V).  (20) 

C  J  C  J 

Test  data  from  REF.  7  (J,  and  C ^  )  for  varying  advance  ratios 
are  shown  on  TABLE  II.  P 

TABLE  II 

Test  (REF. 7)  and  Calculated  data 
varying  loadings 


J 

~ 

T 

c 

CL 

P 

v./v 

J 

v,/v 

A. /A 

J  e 

11.5 

2.87 

5.93 

3.27 

.58 

.5 

3.75 

1.24 

3.48 

2.11 

.64 

.8 

1.18 

.59 

2.12 

1.47 

.74 

1.2 

.37 

.  36 

1.45 

1.15 

.85 

/'Eq?78N 

/e^T? 
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Also  shown  on  TABLE  II  are  flow  and  area  ratios  computed  by  the 
foregoing  method. 

FIG.  5  compares  theory  and  experiment.  At  low  advance  ratios  the 
approximation  for  the  jet  area  permits  closer  agreement  with  the  tests 
than  the  assumption  of  a  noncontracting  jet.  As  the  advance  ratio 
increases  it  is  also  seen  that  the  assumption  of  a  noncontracting  jet 
yields  closer  agreement  with  test  data. 

FIG.  6  is  a  plot  of  the  test  data  J  vs  T  ;  FIG.  7  shows  A. /A  versus 

c  j  ex 

the  thrust  coefficient  For  comparison  with  the  values  on  FIG.  7, 
it  may  be  recalled  that  the  momentum  theory  for  a  free  propeller 
yields  jet  area  ratios  between  1.0  and  0.5  for  zero  and  static  thrusts, 
respectively.  More  work  is  needed  on  the  nature  of  the  interacting 
flow  fields  due  to  the  duct  and  the  exit  flow. 


It  is  interesting  to  note  that  the  magnitude  of  the  lift  coefficients  on 
FIG.  5  results  from  approximately  65%  jet  momentum  component 
and  35%  circulation  effect.  These  values  can  be  obtained  by  compar¬ 
ing  the  vertical  jet  momentum  with  the  total  lift  force  with  power. 


3.  Comparison  with  experimental  data  (REF.  8) 


A  similar  comparison  of  theory  and  experiment  can  be  made  for  the 

data  in  REF.  8.  For  this  source,  the  test  data  notation  are  K.  and 

JL# 

K  ;  the  relation  of  KT  and  K  to  the  coefficients  used  previously  is: 

P  L  p  v  y 

r  1 


L  / 
P/ 


p(wR)2  A 
p(u)R)3  A 


=  CT  (.5\2  /.257Tf) 
■*-'»  P 

=  Pc  (\3/.257Tf) 


(21) 


where  \  =  uR  and  f  =  1  -  (d^/d)2  . 
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0.3 


Angle  of  Attack,  a, degrees 

FIG.  5  LIFT  OF  INCLINED  DUCT  (L/d  =  .5) 
WITH  OPERATING  PROPELLER 


The  geometry  of  the  ducted  system  in  REF.  8  is  as  follows: 

a)  propeller  diameter,  2  feet. 

b)  hub  diameter,  0.43  feet  (d,/d  =  .215),  f  .954. 

n 

c)  exit  diameter,  2  feet. 

d)  duct  length,  6  inches,  duct  no.  P3S. 

e)  L/d  =  .25. 


TABLE  III  below  lists  the  test  data  to  which  Eq.  3  will  be  compared. 

TABLE  III 

Test  data  from  REF.  8  (Fig.  44a, b) 


a 

_X_ 

P  @  ,7R 

kl 

0 

.15 

9° 

.0 

10 

.15 

r 

.0016 

20 

.15 

9° 

.010 

30 

.15 

9° 

.0135 

Inserting  data  back  into  Eqs.  10  and 


K_ 

e 

P 

CT 

_ P 

_E 

c 

L»P 

.0016 

.234 

.355 

0 

.0016 

.234 

.355 

.400 

.0016 

.234 

.355 

.666 

.0016 

.234 

.355^ 

EqMl 

.899 

V./V  =  1.10 
J 


(22) 


Vj/V  =  1.10 

which  seems  to  be  a  low  jet  velocity  value.  The  ideal  momentum 
efficiency  is  .95  based  on  this  jet  velocity  ratio;  however,  in  this 
region  of  jet  ratios  the  influence  of  friction  losses  predominates  and 
low  actual  efficiencies  occur.*  As  before,  the  above  values  are  to 
be  inserted  into  Eq.  3  to  calculate  for  comparison.  FIG.  8  shows 

this  comparison. 


*  REF.  10  The  Flow  Over  Annular  Aerofoils,  by  D.  Kuechmann 

and  J.  Weber.  G.  D.  C.  1 0/ 1 1 33  T,  Ministry  of  Supply, 
British. 
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FIG. 8  shows  a  similar  agreement  between  test  and  theory  as  seen 
on  FIG.  3.  The  duct  L/d  of  the  former  is  .25,  and  is  0.  684  for  the 
latter;  hence,  the  theory,  for  these  isolated  cases  shows  good  agree- 
ment  for  a  fairly  wide  range  of  L/d  values. 


FIG.  8  LIFT  OF  INCLINED  DUCT  (L/d  =  .25) 

4.  Comparison  with  the  experimental  data  (REF.  9) 

The  above  data  allow  an  additional  correlation  between  theory  and 
experiment  for  the  inclined  ducted  propeller.  These  data  are  obtained 
from  a  system  with  the  following  characteristics: 

a)  propeller  diameter,  2.5  feet. 

b)  hub  diameter,  .  625  feet  (d^/d  =  .25). 


-17- 


c)  annular  exit  area,  11%  greater  than  at  rotor. 

d)  duct  length,  1.25  feet. 

e)  L/d  =  0.5. 

REF.  9  presents  data  in  terms  of  the  following  coefficients 

C^,  =  T/pn2d4,  thrust 

0^  =  27tQ/pn2d5,  power 

r]  =  C^J/C  ,  efficiency 

T  p 

CT  =  L/qd  C  ,  lift 
L  ^  s  s 

where  d  and  C  are  the  duct  diameter  and  length  respectively, 
s  s 

TABLE  IV  lists  the  test  advance  ratio  and  thrust  coefficient  at  two 
different  blade  angles.  Corresponding  test  data  on  lift  versus  angle 
of  attack  are  given  on  FIG.  9;  net  efficiency  is  shown  on  FIG.  10. 


TABLE  IV 

test  data  from  REF.  9 


The  blade  angles  of  17°  and  32°  cover  the  maximum  range  for  the 
test  data  at  a  test  speed  of  Mach  0.2. 
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Angle  of  Attack,  a, degrees 


FIG.  9  LIFT  OF  INCLINED  DUCT  (L/d  =  .5) 

WITH  OPERATING  PROPELLER 
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The  procedure  for  calculating  the  lift  to  compare  with  the  test  data  is 
the  same  given  before  with  a  contracting  jet;  namely: 

T  =  pn2  d4CT  =  pAeVe(V.-V)  (23) 

where,  as  before,  V  /V  =  .5(1+V./V);  it  follows  that: 

e  J 

V./V  =  1  +2.448  CT/J2'  (24) 

From  this  point  the  procedure  is  the  same  as  shown  previously. 

FIG.  9  compares  theory  and  test  for  the  lift  coefficient  for  several 
advance  ratios  and  the  two  blade  angles.  The  calculated  lift  coeffi¬ 
cients  are  seen  to  be  in  good  agreement  with  test  data. 

Efficiency  data  on  FIG.  10  are  generally  higher  than  the  other  data 
examined  earlier.  However,  the  efficiencies  on  FJG.  10  should  be 
still  greater  if  only  considered  on  the  estimated  jet  velocities  obtained 
from  Eq.  24.  For  example  at  P  =  17°,  the  calculated  jet  velocities 
vary  from  1.43  to  1.08  over  the  range  of  advance  ratios  and  from  1. 36 
to  1.05  at  P  =  32°,  Ideal  efficiencies  based  on  these  jet  velocities 
approach  20%  higher  efficiencies  than  the  test  data. 

One  explanation  generally  used  to  explain  such  a  difference  is  based 
on  the  fact  that  the  relative  thrust  to  drag  of  the  systen>  is  respon¬ 
sible;  that  is,  lightly  loaded  ducted  systems  can  experience  large  re¬ 
ductions  in  efficiency  from  the  ideal  value. 

D.  Concluding  Notes 

Three  different  set  of  ducted  propeller  lift  data  are  examined  herein 
relative  to  theoretical  predictions.  Agreement  is  generally  good  over 
a  wide  range  of  operating  conditions.  Correlation  at  low  advance  ratio 
need  improvement.  The  behavior  of  the  jet  issuing  from  the  duct  needs 
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clarification  regarding  subsequent  contraction  or  expahsion;  a  calcu¬ 
lation^  of  the  flow  field  of  ducted  propellers  with  appropriate  thrust 
loadings  should  provide  the  needed  results  concerning  the  jet.  Such 
a  study  would  also  be  of  use  for  evaluating  the  inlet  efficiency  with  a 
highly  convergent  inflow  relative  to  the  same  inlet  with  less  inflow 
convergence. 
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II.  POSSIBLE  STATIC  THRUST  IMPROVEMENT 


A.  Relation  of  Thrust  to  Thrust/Horsepower  Ratio 

The  problem  of  improving  static  thrust  can  be  described  by  the 
elementary  theory  which  follows. 

The  static  thrust  of  a  torque-balanced  system*  in  terms  of  radially 
averaged  flow  values  is: 


T  =  pA.V?  (25) 

J  J 

where  A.  and  V.  are  the  jet  cross  section  area  and  velocity  respec- 
J  J 

tively,  both  at  ambient  pressure.  Taking  account  of  the  thrust/horse¬ 
power  ratio  from  Eq.  25  and  8  (with  zero  flow  losses). 


T  1100 

HP  =  V. 

J 


(26) 


Eq.  25  becomes: 

T=  PA .  (^/Hp) 2  (27) 

FIG.  1 1  is  a  plot  of  Eq.  27.  As  seen,  it  is  desirable  to  have  T/HP 
as  large  as  possible  for  a  specified  thrust,  as  well  as  the  smallest 
diameter  consistent  with  good  overall  efficiency. 


Whirl  velocity  component  of  propeller  is  removed  by  a  stator- 
rotor  set. 
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EIG.  11  EFFECT  OF  THRUST /HORSEPOWER 
AND  JET  AREA  ON  STATIC  THRUST 

High  T/HP  values  lead  to  power  economy  as  well  as  low  jet  velocities 
(Eq.  26).  FIG.  11  shows  large  jet  areas  occur  with  increasing  values 
of  T/HP  at  a  fixed  thrust  requirement.  In  principle,  the  derivation 
of  FIG.  11  does  not  require  a  propeller  area  similar  to  the  jet  area; 
the  type  of  duct  and  propeller  design  required  by  theory  to  produce 
a  large  jet  area  with  a  much  smaller  rotor  area  will  be  given  shortly. 

Existing  ducted  propeller  designs  generally  show  values  of  rotor  disk 
areas  similar  to  exit  disk  areas;  also,  the  ultimate  jet  area  usually 
is  ass  »med  about  the  same  as  the  exit  disk  area.  With  such  practice, 
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a  large  jet  area  will  require  a  large  propeller  diameter.  This  point 
will  be  considered  again  later, 

B.  £!low  Values  for  Large  Jet  Area 

A  further  examination  of  the  foregoing  static  thrust  case  is  shown  on 
FIG.  12  where  the  pressure  rise  Ap^  across  the  rotor  disk  is  given 
as  a  function  of  T/HP;  increasing  values  of  T/HP  are  seen  to  require 
decreasing  values  of  Ap^. 

The  derivation  of  FIG.  12  follows  from  Bernoulli  equation 

Po+  Apr  =  Pj  +  (l/2)pV?  (28) 

Eq.  26  furnishes  a  value  of  V.;  since  the  jet  pressure  p.  =  p  , 

j  j  o 

Eqs.  26  and  28  yield: 

Apr  =  (1/2)  p  [  1 100/ (T/HP)] 2  (29) 


0  20  40  60 

T/HP 


FIG.  12  EEFECT  OE  THRUST /HORSEPOWER 
ON  ROTOR  PRESSURE  RISE 
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Since  the  thrust  on  the  propeller  AP^A^  decreases  with  increasing 

T/HP,  (Eq.  29)  it  is  necessary  for  an  increasing  axial  thrust  force 

to  develop  on  the  duct.  This  is  seen  from  the  ideal  ratio  of  duct  force 

F,  to  total  static  thrust  T,  viz:* 
a 

F,  A 

■J-  =  1  -  .5  -j*  (30) 

j 

which  is  obtained  by  subtracting  the  force  on  the  rotor  from  the  net 
thrust.  As  seen,  the  force  on  the  duct  increases  with  increasing  jet 
area. 

To  summarize,  for  a  given  propeller  diameter,  increasing  values  of 
T/HP  require  increasing  jet  areas;  the  pressure  rise  across  the  rotor 
decreases  for  increasing  T/HP-values.  Hence,  the  duct  must  carry 
a  greater  part  of  the  net  thrust.  FIG.  13  illustrates  the  type  of  flow 
needed. 


FIG.  13  EXPANDING  JET  FLOW 


*  The  figure  of  merit  of  a  ducted  propeller  with  an  expanding  jet  is 
2  \lA.iA  1  ;  a  usually  quoted  maximum  value  for  a  ducted  system  is 
2.0.  This  figure  of  merit  is  defined  by  M  =  (T/P)\/T/2  ijA\ 
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Available  test  results  on  a  ducted  propeller  only  consider  ducts  with 
internal  area  expansions;  e.g.,  see  REFS.  11  and  12.  A  rearranged 
form  of  the  data  in  REF.  11  is  shown  on  FIG.  14. 


VT 


Area  Ratio,  A  /A 
e  p 


FIG.  14  EFFECT  OF  INTERNAL  AREA  RATIO  ON  DUCT  FORCE 


Internal  expansion  angles  are  given  on  FIG.  14.  These  data  indicate 
that  flow  separation  from  the  duct  walls  prevented  further  increase  of 
duct  thrust.  Area  ratios  greater  than  the  values  on  FIG.  14  are  pos¬ 
tulated  by  the  flow  type  of  FIG.  13.  The  method  of  REF.  12  seeks  to 
improve  static  thrust  by  large  expanding  internal  exit  areas  with  various 
flow  control  means  to  force  the  flow  to  adhere  to  the  expanding  duct 
walls;  results  do  not  attain  the  promise  of  the  theory  due  to  flow  diffi- 


REF.,  11  Static  Tests  of  a  Shrouded  and  an  Unshrouded  Propeller, 
by  R.J.  Platt,  Jr.  NACA  RM  L7H25.  February  1948. 

REF.  12  Investigation  of  Large  Expansion  Diffusers  For  Minimum 
Rotor  Area  by  J.R.  Duvivier  and  R.  B.  McCallum.  Mass. 
Inst.  Tech.  Aero  and  Struct.  Res.  Lab.  T.R.81-1.  Nov. 

1959. 
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culties  along  the  duct  passageway.  Should  such  flows  be  feasible 
whether  internal  or  external  flow  expansion,  careful  matching  of 
duct  and  propeller  design  will  be  required. 


With  regard  to  the  propeller  design,  FIG.  15  shows  magnitudes  of 
pressure  reduction  p  -  p2  =  Ap2  in  front  of  the  propeller  in  terms 
of  the  pressure  rise  Ap^  across  the  disk  and  varying  area  ratios 
Aj/A^  .  The  pressure  decrease  ahead  of  the  rotor  for  the  expanding 
jet  case  will  always  be  greater  than  the  pressure  rise  across  the 
rotor,  as  seen  from  FIG.  15;  i.e.,  from  Eq.  28,  pQ  =  p^  +  (l/2)pV^ 
and  flow  continuity  between  station  2  and  jet.  Ap  =  p  -  p  . 


FIG.  15  RELATION  BETWEEN  PRESSURE  CHANGES 


Results  on  EIG.  15  follow  directly  from  the  Bernoulli  equation 
applied  to  the  flow  upstream  and  downstream  of  the  rotor. 

The  externally  expanding  flow  of  FIG.  13  may  be  difficult  to  re¬ 
produce  physically.  Jet  boundary  effects  at  the  duct  exit  could, 
significantly  alter  the  ideal  flow  pattern.  However,  it  might  be 
reasoned  that  the  known  occurrence  of  an  expanding  jet  behind  a 
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helicopter  in  the  windmill  brake  state*  is  an  indication  of  possibly 
attaining  a  highly  expanding  jet;  other  known  flows  of  a  related  na¬ 
ture  are  decelerating  inflows  approaching  intakes  and  contracting 
exit  flows  in  an  unbounded  medium.  A  rigorous  analytical  study  of 
this  case  requires  consideration  of  the  duct  shape  and  the  interfer¬ 
ence  effect  of  the  jet  surface  on  the  adjacent  flow  field.  Varying  jet 
shape  problems  generally  are  difficult.  Specific  measurements  re¬ 
lated  to  the  foregoing  discussion  do  not  seem  to  be  available. 

C.  Numerical  Example 

As  an  example  of  the  static  thrust  case  with  expanding  jet,  consider 

a  required  thrust  of  5,000  lbs.  to  be  produced  with  T/HP  =  20. 

FIG.  11  gives  A.  =  695  sq.  ft.  A  12.5  ft.  rotor  diameter  leads  to 
J 

A. /A  =5.7.  Eq.  26  shows  a  jet  velocity  of  55  ft.  / sec .  hence  the 
J  P 

axial  velocity  at  the  rotor  is  314  ft.  /sec. 

Pitch  distribution  for  this  rotor  can  be  shown  to  depend,  with  good 
approximation,  on  the  axial  velocity  at  the  rotor  and  the  blade  rpm. 
FIG.  16  shows  several  radial  variations  of  the  pitch  angle  for  the 
12.5  foot  rotor.  All  the  variations  shown  are  compatible  with  the 
aforementioned  axial  velocity  at  the  disk  and  in  the  jet. 


*REF.  13  See  photographs  in  Report  V.1535,  National  Aero. 
Res.  Inst.  (Amsterdam),  1953. 
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FIG.  16  EFFECT  OF  RPM  ON  DESIGN  VARIATION 
OF  PITCH  ANGLE  ALONG  BLADE  RADIUS 

FIG.  16  is  based  on  a  stator  located  downstream  of  the  rotor;  how¬ 
ever  the  effect  of  an  upstream  stator,  taking  account  of  the  pre¬ 
whirl  velocity,  would  have  only  a  small  effect  on  the  general  trends 
shown  on  FIG.  16. 

The  derivation  of  FIG.  16  is  based  on  the  fact  that 
tan  p  =  TTnd/314 

where  p  is  the  angle  between  the  axial  direction  and  the  relative 
resultant  velocity  approaching  the  blade  profile  .Various  rpm's  are 
selected  and  p  calculated. 

The  pitch  angle  $  is  assumed  to  be  approximated  closely  by  90-P; 

$  is  the  angle  between  the  blade  chord  and  the  plane  of  rotation. 
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Actually  $  =  90-f3+a,  where  a  is  the  angle  of  attack  between  the 
blade  chord  and  the  relative  inflow.  It  can  be  shown  for  the  present 
case  that  a,  generally,  will  be  small  compared  with  9  for  a  blading 
design  based  on  a  potential  or  free  vortex,  type  of  imparted  whirl 
velocity;  see  section  III. 

The  same  power  and  net  thrust  can  be  produced  by  all  the  blade  de¬ 
signs  encompassed  by  the  curves  on  FIG,  16.  Relatively  small 
differences  may  arise  in  the  details  of  the  blade  solidity  and  section 
profile.  Pitch  angle  difference  between  tip  and  root  is  about  20° 
for  all  rpms,  and  tends  to  decrease  slowly  with  increasing  rpm; 
pitch  angles  required,  decrease  with  rotor  rpm  (axial  velocity  is 
the  same  at  all  rpms  on  FIG.  16);  actual  chosen  rpm  would  depend 
on  tip  Mach  number  and  rotational  noise  level;  as  a  guide  to  the  tip 
speed,  Ttnd  =  700  ft/sec  occurs  for  N  -  1069  rpm. 

Pressure  rise  directly  across  the  rotor  for  all  the  above  designs 
will  be  less  than  the  static  pressure  decrease  between  upstream 
ambient  and  rotor  disk  inflow. 

For  comparison  with  the  above  case,  pitch  angles  are  given  on 
FIG.  17  for  a  ducted  system  having  the  same  12.5  foot  rotor,  but 
with  free  jet  areas  decreasing  from  the  695  sq.  ft.  value  of  this 
example.  Shaft  power  for  FIG.  17  is  the  same  as  used  in  FIG.  16.  T 
RPM  for  FIG.  17  is  chosen  at  1400.  The  jet  area  ratio  curve  of 
5.7  on  FIG.  17  is  the  same  as  the  1400  rpm  curve  on  FIG.  16. 

It  is  seen  from  FIG.  17,  that  as  the  free  jet  area  decreases  relative 
to  the  rotor  disk  area,  the  pitch  angle  curves  show  a  lesser  varia¬ 
tion  with  radius  ratio.  For  example,  area  ratio  5.7  at  1400  rpm 
shows  about  a  20°  variation  between  tip  and  hub  compared  to  9°  at 
area  ratio  1.0. 
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1400  rpm 


FIG.  17  EFFECT  OF  FREE  JET /ROTOR  AREA  RATIO  ON 

VARIATION  OF  PITCH  ANGLE  ALONG  BLADE  RADIUS 


The  derivation  of  FIG.  17  is,  as  stated,  based  on  the  same  power 
input  as  the  curves  on  FIG,  16.  This  is  done  by  equating  the  power 
in  terms  of  jet  kinetic  energy;  i.  e. , 


(1/2)  p  (A.V?)b=  (1/2)  p(A.V?)a 


(32) 


where  (b)  denotes  the  values  on  FIG.  17  and  (a)  denotes  the  earlier 

values  on  FIG.  16.  This  equality,  using  (V.)  ~  55  ft/sec,  becomes 

J  a 


(33) 


The  axial  velocity  at  the  rotor  is 


VP  =  <v  Vvj  <34> 


These  values  are  plotted  on  FIG.  18. 
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Knowing  V  ,  the  (3  and  $  angles  can  be  computed  as  denoted  pre¬ 
viously. 

Also  shown  on  FIG.  18  is  the  effect  of  the  free  jet/rotor  area  ratio 
on  thrust/horsepower .  Significant  changes  are  seen  despite  the 
same  rotor  diameter  for  all  calculations.  This  again  implies  the 
duct  design  must  be  compatible  with  these  changes. 
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FIG.  18  EFFECT  OF  FREE  JET/ROTOR  AREA 

ON  THRUST/HORSEPOWER  AND  VELOCITY 
AT  ROTOR 

III.  BLADING  DESIGN  CONSIDERATIONS 

The  previous  section  outlines  requirements  for  investigating  ducted 
propeller  design  possibilities  with  high  thrust/horsepower  ratios. 

In  particular,  the  pressure  rise  across  rotor-stator  set,  force  on 
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duct  and  blade  pitch  distribution  are  given  for  ducted  systems  with 
an  expanding  jet. 

A  further  investigation  of  blading  design  in  terms  of  solidity,  lift 
coefficient,  and/or  the  angle  through  which  each  profile  turns  the 
flow,  is  given  in  this  present  section.  These  results  will  apply  to 
the  blade  design  for  any  type  of  jet:  shape  issuing  from  the  duct.  The 
effect  of  jet  velocity  and  tip  speed  on  blading  design  will  also  be 
shown  in  relation  to  (T/HP)  and  sonic,  tip  speeds  (noise). 

In  order  to  demonstrate  blading  characterics  the  potential  (or  free) 
vortex  type  of  blading  design  will  be  considered. 

The  power  P  absorbed  per  rotor  for  a  free  vortex  design  is  given 
by  the  moment  of  momentum  relation: 

P  =  m  .  TTnd  .  Av  (35) 


where: 

m  is  the  mass  flow  rate  (pAV) 

Av  is  the  tangential  velocity  imparted 
by  the  rotor  tip* 

Eq.  35  can  be  applied  to  flows  with  or  without  losses,  depending  on 
the  manner  for  determining  Av, 

Considering  the  hovering  case,  the  power  without  flow  losses  is 
given  by  the  momentum  relation: 

P  =  i  mV?  (36) 


*  The  product  (d)Av  is  radially  constant  for  free  vortex  blading. 
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which  leads  to 


Av 


,5  V: 

_ I 

ffnd 


(37) 


Eq.  37  is  to  be  combined  with  rhe  often  used  blade  design  relation:* 


crc 


2  |  Av| 
W 

m 


(38) 


where: 

<r 


W 

m 


is  the  solidity  at  a  given  radius 

is  the  lift  coefficient  bas^-d  on  the  mean 
relative  resultant  velocity 

is  the  mean  relative  resultant  velocity 


FIG.  19  illustrates  these  quantities. 

W 


Av 


—4 


TTnd 


FIG.  19  VELOCITY  TRIANGLE  AT  BLADE  TIP 
(without  upstream  stator) 


*  Eq.  38  is  derived  by  equating  the  differential  thrust  on  the 

rotor  expressed  by  the  pressure  rise,  to  the  relation  in  terms 
of  lift  coefficient;  the  pressure  rise  term  is  subsequently  re¬ 
placed  by  the  jet  and  upstream  velocities.  Eq.  38  is  applicable 
with  a  stator  located  upstream  or  downstream.  For  an  upstream 
stator,  combined  only  with  a  rotor,  the  sign  Av  is  taken  as 
negative;  this  accounts  for  the  absolute  sign  shown  in  Eq.  38. 
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It  is  often  more  convenient  for  design  work  to  use  the  stagger  and 
turning  angles  (3  and  9,  respectively,  for  a  chosen  solidity.  How¬ 
ever,  for  evaluating  overall  characteristics  Eq.  38  can  be  very  useful 


The  relative  flow  enters  with  velocity  Wt  and  leaves  the  blade  profile 
at  speed  W2  ,  having  turned  through  an  angle  9  and  imparting  tan¬ 
gential  velocity  Av. 


It  is  convenient  to  use  the  velocity  W  instead  of  W  since  lift  data 

1  l  m 

suited  to  ducted  blading  usually  are  nondimensionalized  with  respect 
to  Wj  .  This  requires  the  change 


W, 


=  c 


m 


U  W 


-r 


(39) 


m 


which  together  with  Eqs.  37  and  38  can  be  placed  in  the  form: 

'2  W, 


crc 


(V./W,)1 
=  (W/W  ) 

'  l  m 


77nd 


(40) 


From  the  geometry  of  FIG.  19,  it  can  be  shown  that 


Wj  sec2^ 

{1T)Z  =  l  tan  -  .25(Aj/A.)‘  cot^  J2  +  1 

J 

and  includes  the  continuity  relation; 


(41) 


V,  =  V.IAVA  )  (42) 

where: 

is  the  annular  area  at  the  rotor  disk. 

A.  is  the  ultimate  jet  area. 
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Eqs.  40  and  41  can  be  placed  into  the  sought  form  which  applies  for 
any  radial  station  having  the  relative  inflow  angle  |3  . 

i 

o-c  =  (V./W)2  {[tan  (3  -  .25  (A  /  A.)2  cot  3]  2  +  l}2  cot  (3  (43) 

J-'1  J  i  1  1  J  1  1 

Eq.  43  contains  only  three  variables,  since  it  can  be  shown  from  the 
geometry  of  FIG.  19  that: 


tan  = 


^-(V/W,)2  (A./A)2 
(A./A1)(V./Wl) 


(44) 


The  velocity  ratio  term  V. /W^  is  of  importance  to  the  ducted  pro- 

J 

peller  for  two  reasons: 


a)  the  jet  velocity  V.  determines  the  (thrust/horsepower) 

J 

ratio  of  the  entire  system.* 

b)  the  relative  resultant  tip  speed  W  governs  the  noise 
level  of  the  ducted  system.** 


Noise  levels  of  a  ducted  system  should  be  much  below  that  of  a  cor¬ 
responding  unducted  propeller.  However,  for  highly  loaded,  high  rpm 
systems,  the  noise  intensity  will  be  disturbing. 


Varying  area  ratios  A  /A.  in  Eq.  43  take  account  of  ducts  with  small 

J 

internal  area  variations  as  well  as  the  highly  diverging  free  jet  ex¬ 
pansion  which  is  predicted  in  section  2. 


* 


It  is  shown  on  page  23  that  T/HP  = 


1100/V. 

J 


It  can  be  shown  that  the  sound  pressure  produced  for  a  given 
harmonic  will  vary  to  a  power  of  the  tip  velocity. 
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FIG.  20  which  follows  is  a  graph  of  Eq.  43. 


FIG.  20  BLADING  CHART  (any  radius) 

The  curve  for  A./Ai  =1.0  approximates  the  blading  design  of  current 
ducted  propellers;  values  somewhat  greater  than  1.0  allow  for  ducts 
with  an  expanding  exit  nozzle.  Area  ratios  of  about  2  and  larger  would 
apply  for  the  free  jet  expansion  described  in  section  II  or  for  expanding 
ducts  having  a  means  for  preventing  flow  separation. 


A  value  of  V7W  of  about  .  1  would  apply  for  V.  and  Wi  in  the  order  of 
100  and  900  ft.  / sec.,  respectively.  Velocity  ratios  near  .  1  lead  to  a 
high  value  of  T/HP  and  a  tip  speed  of  80%  the  sonic  value.  For  V./'WJ 
=  1.  FIG.  20  shows  (rc^  will  be  desirably  low  (<.  10)  over  a  veryJwide 
range  of  A./A^  values.  A  low  <rcL  value  implies  a  lightly  loaded  blade  , 

which  also  tends  to  reduce  the  sound  pressure  generation. 
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As  a  guide  to  the  limit  of  crc.  ,  it  may  be  shown  that  values  for  ac. 
of  2  and  somewhat  greater  is  near  the  upper  limit  for  typical  com¬ 
pressor  blade  profiles  with  high  stagger  angles  (~70°).  However, 
for  detail  design,  the  specific  blade  profile  data  must  be  used. 

While  the  foregoing  analysis  can  also  be  given  in  terms  of  the  turn¬ 
ing  angle  9  shown  on  FIG.  19,  the  presentation  of  FIG.  20  appears 
simpler  for  demonstrating  the  overall  effects  of  velocity,  area  ratio, 
and  velocities  related  to  noise  generation.  It  is  possible  to  relate 
crcj^  to  the  turning  angle  9  and  stagger  angle  |3  at  any  radius.  This 
will  be  done  for  the  stator-rotor  set  shewn  on  FIG.  21. 

Following  the  methods  discussed,  the  relation  for  ox^  becomes 

l 

tan2  ft  -  tan2  (ft  -  9) 

acL  =  - P - t— - -  .  (45) 

1  (tan2  ft  +  l)2  tan  ft 

Eq.  45  is  plotted  on  FIG.  21.  Highly  loaded  profiles,  denoted  by  a 
high  9  and  crc.  -values,  are  seen  to  require  higher  stagger  angles. 
Lightly  loaded  profiles  can  be  accommodated  by  a  much  wider  range 
of  p-values  . 
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FIG. 21  TURNING  ANGLE  CHART 


IV.  ADDITIONAL  DESIGN  CONSIDERATIONS 


A.  Duct 

This  brief  discussion  concerns  practical  considerations  currently 
difficult  to  handle  by  analysis.  Extensive  analytical  work  in  this 
area  is  currently  in  progress  to  simplify  duct  selection. 

It  is  well-known  that  the  flow  field  around  a  duct  without  a  propeller 
is  generally  different  when  the  propeller  is  in  operation.  For  example, 
a  nonpowered  duc^t  having  a  profile  with  the  suction  side  on  the  in¬ 
terior  may  experience  a  greater  or  lesser  suction  with  the  propeller 
in  operation;  such  effects  depend  on  the  design  condition  of  the  en¬ 
closed  rotor.  It  is  also  certainly  to  be  expected  that  the  mutual  in¬ 
teraction  cf  duct  and  rotor  alters  th*  f’ow  field  at  the  enclosed  rotor. 

The  best  choice  for  a  duct  length  is  not  directly  answered  and  can 
depend  on  several  factors;  e.g.,  inflow  acceleration,  pressure  rise 
across  the  rotor,  available  length  required  for  internal  blading  and 
aircraft  stability. 

Ducts  which  are  cambered  to  accommodate  an  accelerated  inflow  may 
experience  inflow  difficulties  at  higher  speeds  where  the  inlet  velocity 
approaches  the  airspeed  and  possibly  a  lower  speed.  For  such  oper¬ 
ation  it  may  be  preferable  to  favor  the  duct  inlet  design  for  the  low 
speed  operation  and  incorporate  a  sufficient  leading  edge  (by  experi¬ 
ence)  radius  to  avoid  intake  flow  separation  at  higher  speeds.  This 
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problem  is  not  always  resolvable  in  a  simple  manner. 

Singularity  techniques  for  studying  ducted  propeller  flows  will  gener¬ 
ally  describe  the  characteristics  of  the  flow  field.  Such  techniques, 
which  generally  follow  the  development  of  REF.  14,  are  cumber¬ 
some  in  application  to  performance  studies  and  it  seems  doubtful  at 
this  time  whether  this  approach  can  be  substantially  improved. 

Current  work  at  ERG  is  seeking  a  rapid  but  highly  accurate  technique 
to  obtain  the  flow  field  for  designing  ducted  systems  utilizing  conven¬ 
tional  techniques  of  combining  mean  camber  lines  with  symmetrical 
profiles . 

B.  Tip  Clearance  Losses 

The  problem  of  selecting  a  tip  clearance  is  one  which  continually 
occurs  with  ducted  propeller  design.  Experience  is  an  important 
guide  for  selecting  a  tip  clearance.  Studies  often  show  the  effect  of 
blade  loading  on  such  losses.  Despite  this,  it  is  interesting  to  note 
the  general  compilation  in  REF.  15  of  tip  losses  from  a  wide  variety 
of  test  data.  A  4.6%  pressure  loss  is  shown  for  each  1%  increase  in 
the  ratio  of  tip  clearance /tip  chord.  The  100%  reference  point  is 
taken  at  a  clearance  ratio  of  1%. 


REF.  14  Fundamentals  of  Annular  Airfoil  Theory  (Nozzles  In  a 

Free  Stream)(  In  German)  by  H.  E.  Dickman.  Ingenieur- 
Archev,  Vol.ll,  1940. 

REF.  15  The  Effect  of  Tip  Clearance  On  the  Peak  Pressure  Rise 

of  Axial- Flow  Fans  and  Compressors  by  L.H.  Smith,  Jr.  , 
AS  ME  Symposium  on  Stall,  1958. 
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C.  Possible  Propulsive  Efficiency  Improvement 


The  externally  expanding  jet  considered  previously  for  V  =  0,  can, 
analogously,  be  shown  for  in-flight  conditions;  this  will  be  shown  for 
the  ideal  propulsive  efficiency  e^;  without  losses  the  familiar  result  in 
terms  of  the  jet  velocity  ratio  is 

e  =  2/[  1  +  V./V]  .  (46) 

P  J 

Eqs.  5  and  7  combine  to  yield 

CT  =  2(A./ A  )(V./V)(V.  / V-l)  (47) 

1  J  e  J  J 


which  yields 


e  = 


p  1.5  +  sli 25  +  ,5Ct/ (A. /A  )' 

x  J  e 


(48) 


For  a  free  propeller  it  can  be  shown  that  A./ A  =  1.0  with  V./V  =  1.0 
r  J  e  J 

and  A. /A  =0.5  with  V  =  0  (static  thrust).  FIG.  22  shows  the  variation 
J  e 

of  ideal  efficiency  with  thrust  coefficient  and  jet  area  ratio.  Efficiency 
is  seen  to  improve  as  the  area  ratio  increases.  Large  effects  are 
shown  as  the  thrust  coefficient  increases.  Design  considerations  to 
explore  this  possibility  are  the  same  as  described  in  part  II  where 
the  expanding  jet  under  static  operation  is  considered. 


-43- 


EASTERN  RESEARCH  GROUP 


m  WALL  SYBEET  NEW  YOU  S.  N.  Y. 


Thrust  Coefficient,  Cr^ 

FIG.  22  EFFECT  OF  EXPANDING  JET  ON  IDEAL  EFFICIENCY 


V.  FLOW  FIELD  STUDIES 


A.  General 

The  analysis  which  will  follow  is  concerned  with  a  jet  producing 
system  located  in  a  duct  and  the  entire  system  located  within  a 
structure  such  as  a  wing  or  fuselage.  FIG.  23  illustrates  the 
general  geometry  of  the  jet  producing  system  and  the  ambient 
air  flow. 

It  is  reasonable  to  anticipate  that  the  inducted  and  ejected  air  flow 
will  produce  interaction  or  interference  flows  over  the  surrounding 
surfaces.  Moreover  the  relative  size  of  the  air  opening  to  the  sur¬ 
rounding  area  can  be  expected  to  significantly  influence  the  magni¬ 
tude  of  the  flow  interference. 

For  example,  with  or  without  an  ambient  stream,  an  inflow  which 
passes  near  to  surfaces  external  to  the  duct  can  alter  the  pressure 
distribution  on  these  surfaces;  a  similar  situation  will  occur  at  the 
jet  exit.  Generally,  unfavorable  effects  can  be  intuitively  reasoned 
for  the  simple  case  of  a  fan-in-wing.  The  issuing  jet  can  be  con¬ 
sidered  as  curved  'tube'.  Upstream  of  this  tube  the  flow  will  be  de¬ 
celerated  similarly  to  a  flow  approaching  a  circular  cylinder;  this 
deceleration  results  in  a  pressure  increase  on  surrounding  upstream 
surfaces.  Lower  pressures  should  occur  around  the  'jet-tube'  be¬ 
cause  the  jet  thickness  and  the  nearby  body  produce  a  venturi  type  of 
effect  to  increase  local  velocities.  The  net  result  of  the  increase 
and  decrease  of  pressures  should  depend  on  the  strength  of  the  vor- 
ticity  related  to  the  jet  velocity  ratio,  but  cannot  be  predicted  with¬ 
out  analysis  or  tests;  however,  the  wing  will  experience  a  probable 
resulting  pitch-up  moment  from  the  pressure  fields,  tending  to  re¬ 
duce  stability  unless  an  additional  compensating  control  is  actuated. 
While  this  isolated  example  can  be  considered  as  a  deleterious  effect 
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of  the  interference  flow,  it  should  equally  be  possible  to  obtain  a 
favorable  effect. 

A  general  mathematical  study  of  the  aforementioned  problem  requires 
a  three-dimensional  theory  of  the  flow  field  around  solid  bodies  with 
interference  flows  produced  by  the  propulsive  system  having  a  finite 
size  jet.  Such  an  approach  will  have  a  wider  application  than  the 
fan-in-wing  case. 

B.  Analysis 

The  analysis  which  follows  treats  the  three-dimensional  problem  of 
a  finite  size  propulsive  jet  interacting  with  surrounding  bodies  of 
finite  thickness.  FIG.  23  illustrates  the  geometry. 


FIG.  23  FLOW  GEOMETRY 
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Formulation  of  the  problem  is  based  on  a  source-disk  flow  issuing 
from  a  nozzle  in  an  annular  disk.  In  principle,  any  surrounding  sur¬ 
face  can  be  chosen.  A  uniform  external  flow  U  is  superimposed  as 
shown.  The  analysis  seeks  the  potential  function  9  of  the  entire  flow 
field  containing  sources,  uniform  flows  and  a  closed  surface;  closed 
surfaces  generally  cause  mathematical  difficulties.  For  spherical 
surfaces,  hydrodynamic al  sphere  theorems  can  be  applied;  this  does 
not  seem  possible  for  the  above  problem.  Hence,  a  different  technique 
is  developed  as  shown. 

Let  $  denote  the  potential  without  the  curved  duct  inlet,  the 

O  u  •  U  • 

contribution  to  $  of  the  source  and  $TT  the  contribution  of  the  uniform 

o  U 

flow.  Hence, 


9  =  9C  +  $ 
o  S.  D.  U 


(49) 


From  the  singularity  distribution  of  a  source  disk,  it  can  be  shown  that: 


-oo 


jm  ,  .  mb  , 

®S.D.,r>z)  1  e 


'slzb  (sb)JJsr)% 

I  V  J 


(50) 


whdre  J's  are  Bessel  functions  and  r  -  \fx*~  +' 


Also,  for  a  uniform  flow: 

^(x,  y,  x)  =  -Ux  (51) 

It  is  known  from  sphere  theorems  in  hydrodynamics  that  for  the 
potential  of  a  system  and  $  =  $q+  the  potential  after  adding  the 
closed  surface  S  to  the  system,  then  $  is  characterized  completely 
by  the  following  three  conditions. 


-47- 


(52) 


a)  $  is  regular  and  harmonic  outside  S 

b)  p2  #  ■♦0  as  pr*’  oo;  (p  =  n/x2  +y2  +  z*  1 ) 

c)  8®  / 8n  =  -8®  /SnonSj  )n  denotes  normal) 

l  o 

For  general  case  of  FIG.  23  the  surface  S  consists  of  the  surfaces 
(see  previous  figure  ): 

(I)  z  =  0,  b  <  r  <  N 

(II)  z  =  a,  b  +c  <  r  <N 

(III)  0  <z  <a,  r  =  N 

(IV)  0  <z  <a  -  c,  r  =  b 

(V)  a-c  <z,  r  <b  +  c,  (r-b-c)2  +  (z-a  +  c)2  =  c2 

Condition  52  c  may,  therefore,  be  written  as  follows: 


8® 

Jz 


(x,  y,  o)  =  (sb)Jo(sr)ds  ;  (b<r  <N) 


8® 


oo 


mb 


-sa. 


9z  (x,  y,  a)  I  e  (sb)Jo(sr)ds;  (b+c<r  <N) 

'O 


a*  .  f*°_  .  . 

-g^(x,y,  a)  =  Ux-~-  /  e  S' Z  'J  (sb)J^  (sN)ds 


'O 

(°<z<a,x2  +  y2  =  N2  ),  where  -^  =  x +  y -g^- 


8®  . 

y,z)  =  Ux-3^ 


C-"" 


!J  (sb)z  ds 

l 


(o<z<a-c,  x2  +  y2  -  b2  ) 


(52c,  I) 


(52cJI) 


(52c,  III) 


(52c, IV) 
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and 


8®  (x,  y,  z)  8«  (x,y,z) 

(r-b-cJ-jJ-  +  (z-a  +c)-j^- 


*»o 


Ux(r-b-c)- (r-b-c)  F  e  8ZJ^  (sb)J^  (sr)ds' 
(z-a 


+  c  )mb  / 

J 

•/ft 


oo 

e  8ZJ  (sb)J  (sr)ds 

i  o 


(52cV) 


[a-c-<z,  r  <b +  c,  (r-b-c)2  +(z-a  +  c)2  =  c2  ] 


Hence,  the  function  4  which  is  sought  may  be  characterized  as  the 
unique  solution  to  (52c,  I)-(52c,  V)  which  is  regular  and  harmonic 
outside  the  surface  S  (condition  52a)  and  such  that: 


lim  p2  ®  (x,  y,  z)  =  0 
p-+  oo 


(53) 


These  are  difficult  equations  to  solve.  The  integrals  appearing  in 
(52c,  I)-(52c,  V)  can  be  evaluated  in  terms  of  elliptic  functions;  how¬ 
ever,  it  is  doubtful  that  such  evaluation  will  simplify  the  problem. 

To  review.the  approach  is  to  let  $  =  be  the  linearized 

o  b.  D.  U 

potential  of  the  system  without  the  closed  surface  S.  _  and®.. 

b  •  D  •  U 

are  both  known  potentials.  If  9  is  the  potential  after  S  is  added,  the 
problem  reduces  to  finding  ®  is  subject  to  conditions  (a), 

(b),  (c)(page  48).  The  problem  can  be  generalized  to  take  angle  of 
attack  a  by  replacing  Eq.  51  by 


®jj(x,  y,  z)  =  -Uxcosa-  uzsina 


(54) 
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The  problem  as  stated  is  to  find  ®  .  Originally  a  solution  was  sought 
by  expanding  Q  in  a  triple  Fourier  series  and  then  determine  the 
coefficients,  presuming  that  a  solution  in  this  form  exists;  whether 
a  solution  exists  or  not,  the  method  involves  cumbersome  calculations. 
For  this  type  of  approach,  a  Fourier  integral  method  now  may  be  more 
appropriate. 

A  more  reasonable  approach  is  to  reduce  the  problem  into  two  simpler 
problems;  namely,  to  find  two  functions  such  that 

(a)  ^  .  is  harmonic  and  regular  outisde  S,  i  =  1/  2 

(b)  p2  .  -*•()  as  pr*  oo,  i  =  1,2 

as  as 

{cl)~5^+~S^  =0  onS 

$  =  $  +  $  is  then  the  required  solution  or  S=S  +  9  +  S  . 

1  11  12  11  12  o 

The  problem  of  finding  involves  only  S  anc  the  uniform  flow  U; 
the  problem  of  finding^  is  circularly  symmetric  about  the  z-axis. 


Potential  $ 

u 


The  problem  is  transformed  to  spherical  coordinates  (r,0,u)  and 
the  equation  of  the  surface  S  is  of  the  form  r  =  f(0).  Its  normal  is 


,1  df  . 

H9)e0 

1+/I 

1  [r  dfl' 


(55) 
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Notation  is  shown  on  the  sketch  below 


In  spherical  coordinates,  if  <|>=ty(r,#,  u) 


x  _  9<t>  ,  1  9<J)  .  1  9<t> 


w  +  Vs\n0  {§7(r2  sinfllf)+l5’(8irf^)+l3(iTS5’^>  ) 


VU(x  cosa  +  z  sina)  =  -Ur  (sin#  cosucos  a  +cos#  sina) 


(56) 


V®  =  - U  {(sin#  cos  ucosa +  cos#  sina) e  +  (cos#  cosw  cosa -sin#  sina  )e  — 
U  r  ; ; 


8®U 

-^=n:  V*TT  =  - 


U 


uiR 


an 

a?' 


-sin  wcos  a} 

[  sin#  cos  (jcos  a  +  cos#  sina-)- 


1  df 

+  —  (cos#  cos  ucosa-sin#  sina  ] 
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Hence  the  boundary  condition  (c2)  for  ftu  on  S  is  (f1  +  df/d 9)  and 


0ft, 

8n 


1 1 


sin  a  (cos  9  -  -p-  sin0  )  (57) 

fl 


+  cos  a  cos  u  (sin  9  +■£*  cos  0)  ]  . 


The  remaining  boundary  conditions  for  ft  appear  previously  as 
(a)  and  (b)  on  page  50. 


A  solution  is  sought  in  the  form  ft  =  F(r,  0)G(u),  +  H(r,0).  Using 
Laplace's  equation  V2ftu  =  0  (in  spherical  coordinates)  it  is  found 
that  sufficient  conditions  on  F,  G,  H  are 


V,2  H  = 

V2  F  = 
2 


d2  G 
dc/ 


+ 


0 

T2  F/  s  in  6 
72  G  =  0 


(58) 


where  7  is  an  arbitrary  constant  and  the  notation  is 

V2  [  <J>(r,0)]  (r2  sin  0  f£)  +  !$-  (sin  9  |^)  .  (59) 

The  boundary  conditions  (57),  in  terms  of  F,  G,  H  on  S  become 

=  U  sin  a  (cos  9  sin0)  (60) 

Sir  fi  flp  fi 

G(-^-  -  ~pr  )  =  U  cos  a  cos  u  (sin0  +  —  cos  9).  (61) 
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Therefore  H  will  not  present  a  difficulty;  set 

00  P  cos^ 


U-UsinaTa  -|+r 

r 


where  the  Pr  are  Legendre  polynomials  and  the  coefficients  an  are 
to  be  determined  such  that  Eq.60is  satisfied.* 


G  and  F  are  still  not  easy  to  determine.  The  expression  which  will 

be  used  is  ^  F  (r0)G  (u)  where  F  and  G  satisfy  the  respective 

y  T  7  yy 

equations  in  Eq.  58.  Hence 


=  b  sin(7i4  +c 

y  y 


cos(7  w) 


To  find  F  the  equationV  2  F=  t2  F/sin0  must  be  solved;  the  family 

y  y  y  +  1 

of  solutions  s\nl9 /y  '  =  F^  are  known.  The  approach  is  to  take 


7  >  1  to  comply  with  the  vanishing  potential  at  least  as  fast  as  1/r2  as 
r  becomes  large.  Hence  the  following  expression  can  be  tried. 


00 

Ucosa  ^ 
m  = 


q  G  (u)  sinn0 /r  n*  * 
n  n 


and  determine  coefficients  qn  to  satisfy  Eq.  61.  After  additional 
solutions  of  the  equation  V2  F=7F/sin0  are  found  it  is  possible  to 
enlarge  upon  this  method  and  obtain  a  wider  variety  of  infinite  series 
available  for  use.  For  example,  there  should  be  a  complimentary 
family  of  solutions  to  the  above  family. 


*  Note  that  for  zero  angles  of  attack,  a  =  0,  it  follows  that  H  =  0. 


Potential  ® 

12 

This  problem  is  axially  symmetric;  cylindrical  coordinates  (R,uz) 

are  to  be  used  as  well  as  the  Stokes  stream  function The  sink 

disk  potential  ®g  ^  is  given  previously  by  Eq.50.  The  equation 

of  the  surface  S  is  previously  written  r=f(0),  now  becomes  ^^fR2-F=  f(arctanR/z) 

which  will  be  written  as  z  =  g  (R).  The  sketch  below  is  a  reminder  of  the 

situation. 


Also  write  dg/dR  =  g*  and 


+1,  z  >  0 

sgn(z)  =  0  ,  z  =  0  =  € 

1 ,  z  >  0 

n  =  (eR  "  g1  ez)  / Vl  +(g1Tn 
v#s  d  [e- Mj^sb^sR.ds) 


e  - 
z 


mb 
'  2 


9®, 


S.  D.  _  r-7  - 

9n  n  •  V^s. D. 


^  J" e’slz[j^  (sbJJ^  (sR)d^e^ 

a??1/ 


oo 

e 


t  (sb)^  (sR)ds  + 

•oo 


/oo 

e  'Z,J  (sb)Jo(sR)ds] 

' 
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oo 

■r-b - r  ("  (  e‘slzljl(flb)[J1(8R)  +  eg1  (sR)]  ds] 

111  +  (g1)2  L  Jo  J 


(65) 


The  boundary  conditions  for  ®  on  S  are  then 


9*  i  0* 
2z  "  8  9R 


■-M 


oo 


-  8  Z 


J,  (sb)[  J,  (sR)  +  eg‘(R)J  (sR)]  ds }  .  (66) 


The  equation  ^  must  satisfy  is 

9z  'R  9z  '  +  9R  'R  8R  ”  u 
noting  that ® constant  as  r  =  ^R*  +  z2  '  -*■  oo . 


(67) 


Since  Eq.  67  is  invariant  under  conformal  mapping,  the  techniques 
of  complex  variables  (two-dimensional  theory)  are  applicable. 

The  fact  that  the  region  S  (more  accurately,  its  intersection  with 
the  R-z  plane)  is  not  connected  causes  some  inconvenience  (it  can¬ 
not  be  mapped  into  a  circle).  More  specifically,  it  seems  clear, 
intuitively,  that  for  the  sketch  above,  a  circulation  is  absent  around 
either  half  of  S;  hence,  there  will  not  be  multiple  valued  functions 
involved  in  the  potential  or  stream  function. 


From  ®  the  potential  ®12  is  obtained  from  the  equations 


9®  9® 

12  1  a*  12 

9R  S  R  9z  ’  9z 


_L  _ 1 

R  9R 


(68) 


The  final  velocity  potential,  as  mentioned,  is 


*  =  *11  +  *12  +  *S.D.  +  *U  • 


(69) 
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C.  Concluding  Notes 


The  analysis  presented  is  believed  to  offer  a  means  to  evaluate 
interaction  and/or  interference  effects  between  the  inflow  and  outflow 
from  ducted  systems  and  surrounding  structures.  It  is  believed 
that  a  digital  computational  program  is  needed  for  rapid  exploration 
of  the  possibilities.  Further,  the  techniques  presented  may  be  used 
to  solve  a  wide  variety  of  related  problems.  A  fan-in-wing  is  one 
such  problem  which  can  be  treated.  For  example,  the  method  could 
be  used  for  the  problem  of  the  reduced  or  increased  pressure  field 
on  the  lower  wing  or  body  surface  which  is  developed  by  the  inter¬ 
action  of  the  fan  exit  jet  and  free  stream  beneath. 
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VI,  APPENDIX:  Derivation  of  Eq.  3 


The  lift  curve  slope  of  a  moderately  inclined  non-powered  duct  is 
given  by  Eq.  2.  With  the  addition  of  a  propeller,  velocities  inside 
and  outside  the  same  duct  will  change.  Moreover,  the  change  in 
vertical  momentum  component  of  the  exit  flow  due  to  power  addi¬ 
tion  must  be  evaluated. 

Singularity  technique  will  be  applied  to  evaluate  the  new  flow  field 
and  lift  coefficient.  The  propeller  may  be  replaced  by  a  disk  of 
uniform  sink  distribution.  It  is  known  that  the  effect  of  such  a 
disk  can  also  be  given  by  a  vortex  distribution  along  the  boundary 
of  the  jet  stream  tube.  The  vortex  distribution  induces  an  axially 
varying,  circumferentially  constant  radial  component  of  velocity 
along  the  duct.  This  radial  component  will  alter  the  duct  shape 
unless  another  vortex  distribution  is  placed  on  the  duct  (cylinder) 
to  induce  an  equal  and  opposite  radial  component. 

Denote  by  a  (x)  along  the  duct  surface  the  inclination  of  a  stream¬ 
line  along  with  the  horizontal  and  represent  the  physical  discon¬ 
tinuity  across  the  duct  surface  by  the  vortex-distribution  y{x) 
along  the  cylinder.*  By  superposition  of  velocities,  the  sum 
of  the  free  stream  velocity  V  and  the  velocity  induced  by  the  vor¬ 
tex  sheet  must  satisfy  the  relation: 

Vr(x) 

V  +  Vx(x)  =  a(x) 


*  The  symbol  7  in  this  Appendix  is  not  related  to  7  utilized  in 
any  preceding  part  of  this  report. 


where  V  and  V  are  radial  and  axial  velocities  induced  by  the 
r  x 

vortices.  For  the  assumption  of  a  slightly  curved  duct  surface 

V>V  and  V  ;  hence  Eq.  70  becomes 
x  r 

V  (x) 

=  a  (x)  (71) 

The  Biot-Savart  law  allows  (x)  to  be  given  by  the  integral 

K7(x)  =  Vr(x)  (72) 

where  y(x)  is  the  vortex  distribution  simulating  the  duct. 

The  Kutta-Joukowski  law  yields  the  side  force  F  on  the  vortex 
sheet;  i.  e. , 

w2 

F=  Jl/2  p  V7  (x)Rdx  (73) 

The  general  approach  is  represented  by  the  foregoing  relations. 

In  the  analysis  which  follows  the  subscript  (u)  denotes  the  flow 
without  power  input;  flow  with  power  will  be  denoted  by  the  same 
symbol  with  subscript  (p). 

The  effect  of  power  input  on  the  lift  will  be  derived  by  first  seek¬ 
ing  the  counter  vortex  distribution  which  maintains  the  same  duct 
shape  defined  by  a(x).  With  power,  it  will  be  assumed  that  the 
flow  inside  the  duct  is  uniform,  axial  and  the  internal  velocity 
V  ^  V  as  the  general  case. 

A  vorticity  distribution  7^  is  introduced  to  simulate  the  effect  of 
the  propeller  on  the  duct.  As  a  result  7^  axial  and  radial  com¬ 
ponents  V^tx)  and  Vr(x)  are  induced  along  the  cylinder. 
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Analogous  to  Eq.  70, 


=  TTvTjO 

x 

where  V  is  an  axial  velocity  along  the  cylinder  which  is  assumed 
to  be  essentially  constant. 

It  is  further  assumed  that  V  is  the  average  of  the  velocities  V 
and  V  outside  and  inside  the  duct;  namely: 

V  =  (1/2)  (V,  +  V)  (75) 

By  virtue  of  email  duct  curvature  and  V»V  (x)  Eq.  74  becomes 

V(x) 

a(x)  =  (1  /2)(VX  +  V)  (76) 

For  a  given  duct  shape,  the  same  boundary  condition  holds  with 
or  without  power;  this  requires  a  (x)  to  be  the  same  with  or 
without  power.  Eq.  76  becomes 


or(x) 


fk7o(x) 

L  v  'J 

u 

Lo/zkv,  +  v)J 

(77) 


Hence  the  relation  between  the  vorticities  with  and  without  power 
is 


7d<*>  *  [  (V,  +  V)/(2V)]  -yo(x)  (78) 

Inserting  7^(x)  into  the  Kutta  Joukowski  relation  (73),  leads  to 
the  relation 


Fd  =  [(Vj  +  V) / (2 V)]  2  Fo 


(79) 
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where  Fq  is  the  side  force  without  the  propeller. 

Eq.  79  does  not  give  the  entire  side  force  on  the  duct;  Fp  is  only 
the  component  due  to  the  vorticity  along  the  duct  surface.  In  addi¬ 
tion.  the  side  force  F.;  is  to  be  added  to  take  account  of  the  vor- 

J 

ticity  along  the  boundary  of  the  jet.  This  force  is  the  difference 
between  the  side  momentum  force  with  and  without  power;  viz: 

F.  =  p(7Td2  /4)(V. V  -V2  )  sin  a  (80) 

J  J  ^ 

for  a  duct  at  angle  a. 


The  total  side  force  F 


P 

F 

P 


with  power  is 

=  F  +  F 
Dr*j 


where  the  lift  coefficient  of  Eq.  3  is 

P 


CT  =  F  / ( 1  / 2 )  pV2  Ld 
P 


and  the  unpowered  lift  coefficient  CT  of  Eq.  2  is 

L,  u 


CT  =  F  /<l/2)pV2  Ld 
■Lit  U  O 


Hence  from  Eqs.  79  through  83,  together  with  the  approximation 

sina=  a,  Eq.  3  for  CT  is  obtained. 

L,  p 


(81) 


(82) 


(83) 
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VII.  SYMBOLS 


A 

B 

C 


I 

m 

n 


P 

Ap 

r 

P 

r 

R 

T 

V 

V. 

J 

Av 

W 


area 

number  of  blades 
blade  chord 
lift  coefficient 

lift  curve  slope 

propeller  diameter 
propulsive  efficiency 

horsepower 
advance  ratio,  V/nd 
loss  coefficient 

duct  length 

mass  flow  rate,  (lb/sec)/g 
rps 

static  pressure 
pressure  rise  across  rotor 

power,  ft  ib/sec 
blade  radius 
d/2 
thrust 

advance  speed 
jet  velocity 

whirl  velocity  imparted  by  rotor 
relative  resultant  velocity 


a  angle  of  attack 

(3  chord  blade  angle  or  angle  of  relative  inflow  with 

axial  direction 

c r  solidity  BC/27Tr 

p  air  density 

X.  L/d  or  uR 

6  turning  angle 

U  27Tn 
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Subscripts 


1  directly  upstream  of  rotor  or  average  velocity  inside 

duct  with  power 

d  duct 

e,  ex  duct  exit 

h  hub 

m  mean  value 

p  powered  or  a  station  in  front  of  propeller 

u  unpowered 


i 
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